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We present preliminary results from a two-pion intensity interferometry analysis from Au+Au 
collisions at ^/s NN =200 GeV measured in the STAR detector at RHIC. The dependence of the 
apparent pion source on multiplicity and transverse momentum are discussed and compared 
with preliminary results from d+Au and p+p collisions at the same beam energy. 

1. Introduction 

Two particle intensity interferometry (HBT) is a useful tool to study the space-time geometry 
of the particle emitting source in heavy ion collisions |l|2j . It also contains dynamical informa- 
tion that can be explored by studying the transverse momentum dependence of the apparent 
source size |3|4| . Extracted parameters in HBT analysis from Au+Au collisions at ^/s NN = 130 
GeV at the Relativistic Heavy Ion Collider (RHIC) did not agree with predictions of hydrody- 
namic models that gave an almost perfect description of the momentum-space structure of the 
emitting source and elliptic flow This "HBT puzzle" could originate from the fact that the 
extracted timescales (emission duration from the R Q /R S ratio and evolution duration from the 
rriT dependence of Ri) are smaller than those predicted by the hydrodynamical model jSJ. 

In this paper we present two-pion correlation systematics as a function of the transverse total 



mass (rriT = ykj- + m 2 , = |(pi + P2)t) and multiplicity in Au+Au collisions at ^/s NN = 
200 GeV produced by RHIC at Brookhaven National Laboratory and measured by the STAR 
detector. We also discuss a fitting procedure in which the strength of the final state Coulomb 
interaction between the two charged pions is taking into account in the fit itself. 



2. Experimental Details 

Experimentally, two-particle correlations are studied by constructing the correlation function 
C*2(q) = A(q)/B(q). Here A(q) is the measured distribution of the momentum difference q = Pi - 
p 2 for pairs of particles from the same event, and B(q) is the corresponding distribution for pairs 
of particles from different events. For this analysis we selected events with a collision vertex 
position within ±25 cm measured from the center of the 4 m long STAR Time Projection 
Chamber (TPC), and we mixed events only if their longitudinal primary vertex positions were 
no farther apart than 5 cm. We divided our sample into six centrality bins, where the centrality 
was characterized according to the measured multiplicity of charged particles at midrapidity. 
The six centrality bins correspond to 0-5% (most central), 5-10%, 10-20%, 20-30%, 30-50% and 
50-80% (most peripheral) of the total hadronic cross section. Charged pions were identified by 
correlating their specific ionization in the gas of the TPC with their measured momentum (> . 

The effects of track-splitting (reconstruction of a single track as two tracks) and track-merging 
(two tracks with similar momenta reconstructed as a single track) were eliminated as described 
in [7]. A new procedure to take into account final state Coulomb interaction is described in the 
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next section. 

The effect of the single-particle momentum resolution (<5p/p ~ 1% for pions) induces system- 
atic underestimation of the HBT parameters. Using an iterative procedure we corrected our 
correlation functions for finite resolution effects. The correction due to the uncertainty on the 
removal of the artificial reduction of the HBT parameters associated with the anti-merging cut 
has been calculated in |Sj and is included as systematic error. 

3. Fitting procedure 

The three-dimensional correlation functions were generated. The relative momentum was 
measured in the longitudinal co-moving system (LCMS) frame, and decomposed according to the 
Pratt-Bertsch |9|10] "out-side-long" parametrization. There is a Coulomb interaction between 
emitted particles that needs to be taken into account in order to isolate the Bose-Einstein 
interaction. This Coulomb interaction, repulsive for like-sign particles, causes a reduction on the 
number of real pairs at low q, reducing the correlation function. In our previous analysis [Till j 
as well as in previous experiments, this was corrected by applying a pair Coulomb correction to 
each pair in the background |H] corresponding to a spherical Gaussian source of a given radius; 
we call this standard procedure. The correlation function was then fit with the functional form: 
C(q a ,q s ,qi) = £(q)xl^i(g im ,) = 1 + A exp(-Rl q z -R 2 s q 2 s - Rf qf), where K coul {q inv ) is the square 
of the Coulomb wave-function. However, this procedure overcorrects the correlation function 
since all background pairs are corrected, including those that are not formed by primary pions. 

We have implemented a new procedure, first suggested by Bowler ^2] and Sinyukov |13| 
and recently used by the CERES collaboration [14] . in which the strength of the Coulomb 
interaction is taken into account in the fit itself and only pairs with Bose-Einstein interaction 
are considered to Coulomb interact; we call this Bowler- Sinyukov procedure. The fit in this 
case is: C(q , q s , qi) = = (1-A) + A K cou i(q inv )(l + exp(-R 2 q z . #2 q 2 _ R 2 w h e re 
Kcoui(qinv) is the same as above. 

In Figure^ the measured 7r + 7r~ correlation function is compared to several calculations. Lines 
indicate the standard {K cou i{Qi nv )) and Bowler- Sinyukov (1 + (A — l)K cou i(Qi nv )) Coulomb 
functions; in the latter, A was extracted from the fit to the 3D like-sign correlation function. 
Clearly, the Bowler- Sinyukov function better reproduces the data. Further improvement is 
observed when the strong interaction (negligible for like-sign pion correlations) is included ^S] 
into the ir + ir~ final state interactions. 

When we use this procedure in our 3D analysis we observe an increase in R Q of 10-15%. The 
values of R s and Ri do not depend significantly on the Coulomb procedure. Consequently the 
increase in R Q /R S is not enough to solve the HBT puzzle. 

4. HBT parameters versus centrality and transverse momentum 

Figure |2] shows the dependence of the source parameters for pions at six centrality bins 
from Au+Au collisions as well as from p+p and d+Au collisions at same beam energy. The 
three radii increase with increasing centrality and Ri varies similar to R Q and R s ; for R Q and 
R s this increase may be attributed to the geometrical overlap of the two nuclei. The extracted 
radii rapidly decrease as a function of uit, which is an indication of transverse flow [T2|. In 
order to extract information about the dependence on centrality, we fit the tut dependence 
of each radius and each centrality to a power-law function: i^(my) = Rio ■ m^ a (solid lines 
in Figure I2J). Figure 01 shows the dependence of a on the number of participants; for Au+Au, 
a is constant for R[ as a function of number of participants and decreases with the number of 
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Figure 1. Experimental (triangles) and the- 
oretical |T3] (stars) ID ir + ir~ correlation 
functions compared with standard (discon- 
tinuous line) and Bowler- Sinyukov (contin- 
uous line) functions. 



Figure 2. HBT radii for 6 different cen- 
trality from Au+Au collisions, and from 
p+p and d+Au collisions. The lines indi- 
cate power-law fits to each parameter and 
centrality. 



participants for R Q and R s for the most peripheral bins indicating a decrease of transverse flow 
for these collisions. R a /R s ~ 1 which indicates a short emission duration in a blast wave fit 17 . 
Assuming boost-invariant longitudinal flow we can extract an evolution time-scale by using 

a simple fit ^7]: Ri =< tj > ^t^ ^t/tO wnere T is the freeze-out temperature and K\ 

and K2 are the modified Bessel functions of order 1 and 2. For T extracted from fits to pion, 
kaon, and proton transverse momentum spectra (90 MeV for most central and 120 MeV for 
most peripheral collisions) ^S] we get < tf a >~ 9 fm/c for central events and < tf >~ 6 fm/c 
for peripheral events. Hence, the evolution time, in addition to the emission duration, is quite 
short. 

For a transverse expanding, longitudinally boost-invariant source, and assuming a Gaussian 
transverse density profile, we can extract information about its radius, R ge om, by fitting the 

I R? 

rriT dependence of R s to: R s (rriT) = J 1 1 ^ ^t°™m T ^ where T is again the freeze-out temperature 

and rjf is the surface transverse rapidity |19|. For T and rjj consistent with spectra we see an 
increase on this radius from ~5 fm for the most peripheral case to ~13 fm for the most central 
one as shown in Figure |1I We also observe a smooth transition from p+p (N par tici P ants = 2) and 
d+Au (iV participants = 8.3) to Au+Au collisions. 

5. Conclusion 

We have presented identical pion interferometry results for Au+Au collisions at ^/s NN =200 
GeV. With respect to multiplicity and uit dependencies, pion HBT radii are very similar to 
results reported at ^/s NN =130 GeV. HBT radii and geometrical radius increase with increasing 
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Figure 3. Extracted a parameter from the 
power-law fits to the HBT radii (lines in Fig- 
ure [2j for p+p, d+Au and 6 different central- 
izes in Au+Au collisions 



Figure 4. Extracted R ge om radius for a trans- 
verse gaussian density profile ^5] for p+p, 
d+Au and 6 different centralities in Au+Au 
collisions 



centrality. HBT radii decrease with and we observe a stronger flow for the most central 
collisions. Our results indicate that both the evolution timescale (as measured by the 
dependence of R{) and the emission duration (probed by comparing R a to R s ) are surprisingly 
fast. The Bowler-Sinyukov Coulomb procedure does not solve the "HBT puzzle" although 
increases the ratio R /R s by 10-15%. 
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